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^ , Abstract 

(— I . We prove upper bounds on the ground state energies of the one- and two-component 

d • charged Bose gases. The upper bound for the one-component gas agrees with the high 

density asymptotic formula proposed by L. Foldy in 1961. The upper bound for the two- 
component gas agrees in the large particle number limit with the asymptotic formula 
conjectured by F. Dyson in 1967. Matching asymptotic lower bounds for these systems 
were proved in references I IQJ and jJJ. The formulas of Foldy and Dyson which are 
based on Bogolubov's pairing theory have thus been validated. 



1 Introduction and main results 

In 1961 L. Foldy [7] used Bogolubov's 1947 pairing theory [^ for Bose systems to give a 
heuristic calculation of the ground state energy of a one-component charged Bose gas in the 
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high density hmit. The one-component Bose gas is a system of Bose particles all of the same 
charge moving in the presence of a fixed uniform background of the opposite charge. 

In 1967 F. Dyson [H] considered the two-component Bose gas with two species of bosons 
with opposite charges. Motivated by Foldy's calculation Dyson was able to prove a rigorous 
upper bound on the ground state energy. A famous consequence of Dyson's upper bound 
is that charged bosonic matter is not stable, the ground state energy is super-linear in the 
number of particles. Dyson, moreover, conjectured an exact asymptotic form of the ground 
state energy in the limit of a large number of particles. 

In [TUl it was proved that Foldy's calculation is indeed correct as a leading asymptotic 
lower bound for the ground state energy of the one-component charged Bose gas in the high 
density limit. 

In pT| it was similarly proved that Dyson's conjectured expression is correct as an asymp- 
totic lower bound for the ground state energy of the two-component charged Bose gas in the 
limit of a large number of particles. 

The aim of the present paper is to prove the corresponding upper bounds thus validating 
both Foldy's one-component and Dyson's two-component formulas. 

It should be mentioned that Foldy's calculation may be viewed as a trial state calculation 
and may thus be turned into a rigorous upper bound. Foldy, however, uses periodic boundary 
conditions, and a periodic version of the Coulomb potential. It is not known whether this 
formulation has the same thermodynamic limit as the formulation given below. 

The one-component Bose gas is a system of A^ particles all of the same charge -|-1, say, 
constrained to a box A = [0,L]^ C M^, in which there is a uniform background charge of 
density p. 

The Hamiltonian for the one-component charged Bose gas is thus 

N 

^iJ^ = E(-l^^-^(^^))+ E \^^-^^-'+c (1) 



where 



^2 



V{x) = p j \x^-y\ ^ dy, C = — \x - y\ ^ dx dy. 

J A ^ JJaxK 

We use Dirichlet boundary conditions. 

It is known from the work of Lieb and Narnhofer |9J that the ground state energy E^^\N) 
of Hpj has a thermodynamic limit if we restrict to a neutral system 

EW{N) 

<p) = 1^-^ — ^^- 



It is however also shown in jU] that one will get the same thermodynamic energy by mini- 
mizing over all particle numbers, i.e., 

e(p) = lim inf — — . 

L^oo N L^ 

THEOREM 1.1 (Foldy's formula). The ground state energy e(p) of the one-component 
charged Bose gas satisfies the asymptotics 



-5/4 

where 



lim p-'/'e{p) = -Jo, (2) 

p~*oo 



The two component Bose gas is described by the Hamiltonian 



H 
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(2) _ Y- ^ A , ^:^ ^i^j 



E-^^.+ E 



N / J <-) — ' ' / J I 

^ Ju 7 ^ •J I 

j=l l<i<j<Af -" 

acting on the Hilbert space L^(R^ x {1, —1}), where the variable (xj, Cj) G M'^ x {1, —1} gives 
the position and charge of particle i. 

The word tvjo component refers to the fact that the charge of each particle can be either 
positive or negative. Thus the gas has a positive and a negative component. One would 
not normally consider the charges as variables, but rather fix them to have given values. If 
we did that, the Hamiltonian would not be fully symmetric in all A^ variables, but only in 
the variables for the positively charged particles and negatively charged particles separately. 
Clearly, the charge variables commute with the Hamiltonian and the bottom of the spectrum 
(the ground state energy) E'^'^^N) of Hj^ will therefore be achieved for a fixed combination 
of charges (rather than a superposition). 

THEOREM 1.2 (Dyson's formula). The ground state energy E^'^^N) of the two-component 
charged Bose gas satisfies the asymptotics 

lim N-'^/^E'^^\N) = -A 

where A is the positive constant determined by the variational principle 



A = infU / |V<l>P-/n / <l>^/2 



< $, / $2 = 1 I (4) 



with Iq again given by ^. 



In |n] Dyson proves that E^'^\N) < —CN"^^^, but with a constant different from A. He 
conjectures that the correct value is given as above. That the exponent 7/5 is, indeed, 
correct was first proved in 1988 by Conlon, Lieb, and Yau in |2|, where they show a lower 
bound —CN"^^^, but still not with the correct constant. They also proved that 5/4 is the 
correct exponent in Foldy's formula. The asymptotic lower bounds in Theorems II. II and II .21 
were proved in ^U] and ^I] respectively. The main results of the following paper are the 
asymptotic upper bounds. 

In Sect. 121 we give a general construction of bosonic trial states on the bosonic Fock space 
over a general Hilbert space. The trial states will be build from coherent states and squeezed 
states. The trial states are essentially the ones dictated by Bogolubov theory. These trial 
states are the bosonic equivalent of the fermionic states in Hartree-Fock theory or rather to 
their extension including the Bardeen-Cooper-Schrieffer states (see pH). 

In the same way as fermionic systems may be approximated by the semi-classical Thomas- 
Fermi theory we will also use a semi-classical type approximation to the Bogolubov trial 
states. 

In Sect. El we use the general trial state method to give an upper bound on the ground 
state energy for the two-component gas, but in a grand canonical setting where we do not 
fix the total number of particles. 

In Sect 13. II we show how to get an upper bound for fixed particle number and thus prove 
Theorem 11.21 

In Sect. Elwe use the general trial state method to give an upper bound on the ground 
state energy for the one-component gas and prove Theorem 11.11 

A key ingredient in the proofs is a semiclassical construction where we represent operators 
as phase-space integrals with coherent states symbols and use the Berezin-Lieb inequalities. 
We need an operator version of the inequality. This is discussed in Appendix ^ 

Acknowledgment: I would like to thank Elliott Lieb, Kumar Raman, and Robert 
Seiringer for valuable discussions. 

2 The abstract trial state construction 

Our goal in this section is to construct trial states on the bosonic Fock space JF = jF(7ii) = 
©^=0^^^' ^^^^ some Hilbert Space Hi, i.e., Hn = (S^sym'^i ^^^ '^o = C- We will be using 
the language of bosonic creation and annihilation operators as a convenient tool for the book 
keeping. We denote by |0) the vacuum vector in JF. If T is an operator on ?ii and W is an 
operator ?-^i ® Tii, which is symmetric under interchange of the tensor factors, we may lift 



(also referred to as second quantize) these operators to T as 

OD N oo 

QY,T. and ® Y. »'«■ 

N=l i=l N=2 l<i<j<N 

Here Tj refers to the operator T acting on the ith factor in the tensor product and Wij refers 
to W acting on the ith and jth factors. If u^, a = 1, ... is an orthonormal basis for ?ii we 
can express these operators using creation and annihilation operators as 

oo TV 

^ ^ Ti = ^(wa, Tup)a{uaya{up) (5) 

JV=1 1=1 a,l3 

and 

oo _. 

® ^ Wjj = ^ y^ (mq U/3, WUf, ® Uu)a{ua)* a{upy a{uu)a{u^) . (6) 

Ar=0 l<i<j<N ap^iv 

Of special interest is the number operator (the second quantization of the identity) 

oo 
N=0 

If G Til is a not necessarily normalized vector we define the corresponding coherent state 
as the normalized vector in Fock space 

|0)c = exp(-||0||V2 + a(0)*)|O) 



oo 



= 5:eHI*/.W)T|0). ,7) 

n=0 "• 

and for a normalized ip G Tii we define the squeezed state depending on A G C with |A| < 1 

|A;^)s = (l-|AnVSxp(-(A/2)a(V^)*a(V^)*)|0) 

= (l-|Ap)^/^f;^^^(a(V^)*)^lO). (8) 

n=0 

It is straightforward to check that these states are normalized. Up to an overall phase \(f))c 
and |A;-?/')s are characterized by 

(a(0)-||0f)|0)c = O and (a(V^) + Aa(V^)*)|A; V^)g = 0. (9) 

We immediately see that 

c(0|(a(0)*)'"a(0)'=|0)c = ||0f ('"+'=). (10) 



For the squeezed state we get 



2n + : 
(n + k) 



l-W')"'i;^l^(2"-^ + l)(2"-^ + 2)-(2") 



n=0 

2n/ \ /r,\k 



x{n + k){n + k-l)---{n + l)(|A|/2)^"(-A/2)* 

- a-|A|V^|Art-A)'-i^(|A|-'^)\l-|An-"^ (U) 

Moreover, the expectation in the state \X;iIj)s oi a product of an odd number of the operators 
a{tp)* or a{ilj) vanishes. 

For the expectation of the particle number we find 

c(0|a(0/||</)||)*a(</'/ll0ll)l0)c = Ur and s(A; ^|a(^)*a(^)|A; ^)s = -^L. 

1 — |Ap 

We point out that the variation in the particle number is very different in the coherent state 
and in the squeezed state 

c(0|(a(0/||0||)M0/ll0ll))V)c-c(0|a(0/||0||)M0/ll0ll)l0)c = Uf (12) 

s(A;7A|(a(^ra(^))'|A;^)s-s(A;^|a(^)*a(^)|A;^)i = ^'^|' (13) 

Thus in the coherent state the standard deviation of the particle number is the square root of 
the expectation itself, whereas for the squeezed state the standard deviation of the particle 
number is, in fact, greater than the expectation itself. For this reason the squeezed states 
are not appropriate for describing Bose condensates with a macroscopic and sharply defined 
occupation number in a specific one-particle state. To describe condensates we will use 
coherent states. 

We will here define a variational principle corresponding to the Bogolubov theory of Bose 
gases. We shall do this by characterizing the set of variational trial states. 

The Bogolubov variational theory is very similar to the Hartree-Fock theory for Fermi 
gases. More precisely, it is similar to the generalized Hartree-Fock theory which includes the 
Bardeen-Cooper-Schrieffer (BCS) trial states. In generalized Hartree-Fock theory (see ^J) 
the class of trial states is defined to be the quasi-free states on a fermionic Fock space. For 
the ground state (zero temperature) theory we may restrict to pure quasi-free states. 

To describe the variational states of Bogolubov theory we we will again start from (nor- 
mahzed) quasi- free pure states. Such a state may be characterized as follows. If \E' G J-'(TCi) 
is a normalized quasi- free pure states there exists an orthonormal family ipi, . . . of Tii and a 
sequence of numbers < Ai, . . . < 1 with Yl'^=i -^q < °^ ^^^^ ^^^t 



*=n 



0=1 



1 - A^)4 exp ( — ^a(z/'Q)*a(z/'o 



|0) (14) 



A straightforward but lengthy calculation from (jllj) shows that the quasi-free state satisfies 

(^, 4444^) = (^, 44^) (^, 44^) + (^, 44^) (^, 44"^) 

+(^,44^)(^,44^) (15) 

and from the definition of the state we have for all integers m > 1 

(vi/,4...4^„,vI/) = 0. (16) 

In (jl5|) and (|16|). a^-, j = 1, 2 . . . refer to any creation or annihilation operators. The relation 
(jl5|) is the case m = 2 oi the more general rule 

(^, 4 ■ ■ ■ «L^) = E (^' «l(i)«l(2)^) ■ ■ ■ (^' «l(2™-i)«l(2^)^), (17) 

7reP2m 

where P2m is the set of pairing permutations 

P2m = {vr G ^2^ I 7r(2j - 1) < 7r(2j + 1), j = 1, . . . ,m - 1 

7r(2j-l)<7r(2j), J = l,...,m}. (18) 

We shall here use this only in the case ()15p when m = 2. 

The one-particle density matrix of the quasi-free state \1' is the operator 71 defined on 
the one-body space Hi by {g,lif)ni = (^) 'j^(/)*'j^(fl')^)jp' where f,g E Hi. From (fTT|) 

00 .2 

Note, in particular, that the one-particle density matrix is a positive semi-definite trace class 

operator with 

00 ,2 

Tr7i = {^,Af^) = J2 T^ry < ^• 

a=l " 

Connected to the quasi-free pure state \E' we also have the symmetric bilinear form ^1 on Hi 
given by ^i{f,g) = {'^,a{f)*a{gy'^)yr. We find, again from ((TTI), that 

00 , 

a=l " 

We may identify ^1 with a linear map ^i : Hi ^ HI, from the one-body space Hi to its 
dual space HI- We then have the relations 

^1*6 = 71(71 + 1) 671 = 716, (21) 



where we have also identified 71 in the natural way with a map from 7i* to itself. If we 
introduce the operator T : Hi Q) HI ^ Hi Q) HI defined using matrix notation as 




we may rewrite the condition ()21|) as 




r. 



We may refer to an operator satisfying this condition as a symplectic projection. In the 
fermionic case the corresponding operator is simply a projection. Note that the operator F 
may also be described by 

(l/i) © (^il, ri/2) © {92\)n,^ni = (^' W^r + <92)){a{h) + a((7i)*)^)^(w,) , 

where we have used the Dirac bra and ket notation to denote elements of Hi and H\ respec- 
tively. 

Given a positive definite trace class operator 71 and a symmetric bilinear form ^1 satisfying 
fl21|) we may find a unique quasi-free pure state \1/ such that 71 is the corresponding one- 
particle density matrix and ^1 the corresponding bilinear form. To see this one simply has 
to show that there exists an orthonormal family ipi, . . . and a sequence of positive numbers 
Ai, . . . such that (J19j) and (j2Up hold. This is a fairly simple exercise in linear algebra. 

The choice of ^1 is equivalent to a particular choice of eigenbasis for 71. If 71 has real 
eigenfunctions (in some representation) there is a particular ^1 corresponding to this choice 
of basis. We shall use this in our construction of states in the next sections. 

Consider as an example 71 being a real translation invariant operator on the Hilbert space 
L^(]R"/27rZ") of square integrable functions on the torus. The real eigenfunctions come in 
degenerate pairs of the form cos{px) and sin(px), p G Z". The associated quasi-free state 
will in the exponent have terms of the form 

a(cos(px))*a(cos(px))* + a(sin(px))*a(sin(px))* = a(e'P'')*a(e"*^'')*. 

This corresponds to a pairing of states with opposite momenta, as is the usual case in the 
Bogolubov pair theory. 

The Bogolubov variational states are not just quasi-free states as defined above. In 
fact, quasi-free states being build out of squeezed states are not well suited for describing 
condensates (see the discussion after (fT^ and (jT^. We introduce condensates by appropriate 
unitary transformations of quasi-free states as we shall now describe. 

8 



Given cj) G Tii we have a unitary map U^ on the Fock space Til-Li) which satisfies 

U;a{f)U^ = a{f) + {f,ct>). 

This unitary is unique up to an overall complex phase, which we may fix by noting that we 
can add the requirement that the unitary maps the vacuum state to a a coherent state 

UM = \4>)c- 

From the first identity in Q it is clear that f/^ satisfies this up to a phase. 

The Bogolubov variational states are constructed from a quasi-free state \1/ and a vector 
G Til as \l/0 = U(j,^ ■ From the above discussion we see that a Bogolubov state may be 
described as follows. 

DEFINITION 2.1 (Bogolubov variational states). A Bogolubov state on the bosonic 
Fock space J-" (Hi) is given by 



^<^m,Ci = n (1 - A^)iexp [-^(ai^^y - (0,^„))(a(^„)* - {cP,^^.))] 



:22i 



a=l 

where (f) G Tii and ipi, ip2 ■ ■ ■ is an orthonormal family in Tii and < Ai, A2, . . . < 1 satisfy 
S^i -^a = 1- We call (p the condensate vector and ipi,ip2 ■ ■ ■ the pair states. 

There is a one-to-one correspondence between Bogolubov states and triples (0, 7i,,^i) 
consisting of a vector G TYi a positive trace class operator 71 on Hi and a bilinear form ^1 
on Hi X Hi satisfying pTjl . The correspondence is given by (fTn|l and pUjl . 

We find for the one-particle density matrix of the Bogolubov state \l/^,^i,^i that 

(^0,7i,a' «(«)*«(^)^</',7i,6)^(Hi) = (^0,71,6' («(«)* + (</'' ^))(a(^) + (^' </'))^0,7i,Ci)^(7^i) 

= (t;,7iw) + (t;,0)(0,w) (23) 

and likewise for the two-particle density matrix using (|15|) 

(^</',7i,5i' «(^i)*«(«2)*a(w2)a(^^i)^<?i,7i,5i)^(?^i) = (^'i' </')(^2, 0)(0, Wi)(0, ^^2) 

+6 («1, ^2) (^1, 0) (^2, 0) + ^1 (Vl, V2) (0, Ml) (0, M2) 
+ (t;2,7lMl)(t;i,0)(0,M2) + {vi,-fiU2){v2,(t>){(t),Ui) 
+ iv2,'JlU2){Vl,(p){<P,Ui) + (t^l,7l^il) (^2,0) (0,^2) 



+ (^^l,7lMl)(t^2,7lM2) + (fl,7lM2)(f2,7lMl) +^lK,f2)6(Ml,M2)- (24) 

The above trial states are motivated by the Bogolubov approximation for Bose condensed 
systems. The states represents the condensate, whereas the states -00, « = 1, • • • represent 
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the pair states. A key ingredient in the Bogolubov approximation is the c-number substitu- 
tion, i.e., the replacement of the operator a(0) by the number ||0|p. This replacement will 
give the correct value for expectations of normal ordered products in the Bogolubov states 
if we have the additional assumption that 7i0 = (see (fTUj) . In Section El we will choose a 
Bogolubov state satisfying this assumption, but in Section HI the Bogolubov state that we 
choose will not satisfy the assumption. 

It is not the aim here to study the general properties of the Bogolubov variational prob- 
lem, i.e., the minimization of the expectation of many-body Hamiltonians restricted to Bo- 
golubov states. We will instead proceed to the specific examples of the one-component 
and two-component charged Bose gas. Here we shall not characterize the exact Bogolubov 
minimizer, but instead give the semiclassical approximations to these states which give the 
leading order asymptotics in Theorems 11.11 and 11.21 

The Hamiltonians that we are interested in are particle number conserving, i.e., commute 
with particle number and the reader may wonder why we do not define a class of particle 
conserving, i.e., canonical trial states rather than the grand canonical states above. As in the 
fermionic BCS theory it is very complicated to write a canonical trial state. The calculations 
are greatly simplified in the grand canonical setting. Simple minded trial states with a fixed 
number of particles in the condensate will not give the correct approximation, since the 
important virtual pair creation will be lost. 

3 The two-component charged Bose gas 

We consider the two component Bose gas described by the Hamiltonian 

N 



h'"-®hs\ j?^='=e4^'+ e 



CiCj 



— / V 

N=0 i=l i<i<j<N ' "" 

acting on the Fock space J^{L'^{M.^ x {1,-1}), where the variable (xj,ej) G M'^ x {1,-1} 
gives the position and charge of particle i. 

Our goal here is first to construct a grand canonical normalized trial function 

^ eJ^{L\R^x {1,-1}) 

with particle numbers concentrated sharply around the average value (A/") = (\E',A/'\E') and 
such that 

(i/(2)) = (^, i7(2)^) < -A{J\fy/'' + o{{U)'''^) (25) 
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for large (A/"). We have denoted the expectation in the state \1/ by (A) = (\1/, A"^). From this 
the proof of Dyson's formula Theorem 11.21 (i.e., the fact that we can achieve this estimate 
with a trial function of fixed particle number) will follow fairly easily (see Section IXT|) . 

To construct the trial state \E' we use the method from the previous section. We begin 
with a normalized minimizer <l> for the variational problem (@)) . Using spherically symmetric 
decreasing rearrangements it is not difficult to see that a minimizer exists and that it may be 
chosen positive and spherically symmetric decreasing. Moreover, from the Euler-Lagrange 
equation it is exponentially decreasing and smooth. It is, however, not essential that we can 
find an exact minimizer with these properties. As we shall see, we could as well have chosen 
an approximate minimizer, which is smooth and compactly supported. 

Let n > and define the normalized function 

0o(a;)=n3/io$(ni/5x). (26) 

We define a normalized state \E'„ € JF as in (j^^ with the condensate vector on L^(R^ x 
{-1, 1}) given by 

71 

0(x,e) = J-(f)o{x) 
and the operator 71 on L^(]R^ x { — 1,1}) defined by the integral kernel 

7i(a;,e;,|/,e') = --f{x,y)ee', 

where 7 is a positive semi-definite trace class operator having real eigenfunctions. We shall 
make an explicit choice for 7 below fseelS^. We write the spectral decomposition of 7 as 

0=1 '^ 

where ^/'q, a = 1, ... is a real orthonormal basis and < Aq, < 1 for a = 1, . . .. Observe that 
on the space L^(]R^ x {1, —1}) we have ||0|p = n and 7i0 = 0. Denoting 

^a±{x,e) =tjJa{x)5±i^e, a = 1, . . . 

we may write the trial state \l'n as 

^n = n(l - AD'/'exp ( -^ + a*(0) - E E ^^^'<ea:A |0) , (28) 

a=l \ e,e'=± a=l / 

where a* g = a{ipae)*, for a = 1, . . .. 

As discussed in the previous section choosing n and any 7 with real eigenfunctions 
uniquely specifies a state \E'„ of the form above (possible degenerate eigenvalues will not 
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cause ambiguities). Instead of specifying the individual eigenfunctions ip^ and parameters 
Aq,, a = 1, ... we will simply choose the operator 7. 

The state \l/„ should be compared to Dyson's trial state in [0]. The main difference is 
that whereas we use a coherent state construction for the condensate, Dyson used squeezed 
states for this as well. Put differently, Dyson's trial state corresponds to an exponential of a 
purely quadratic expression in creation operators without any linear terms. As we explained 
in the previous section the consequence of using the linear term in the exponent is that 
the variation in the number of particles occupying the state 0o is much smaller than for a 
quadratic term. 

From (j23j) we find for the expected number of particles in the state ^n 

(■^) = ( E E <e^-e ) = ri + Tr7. (29) 

and for the kinetic energy expectation 

00 N 



N=0 i=l 



AT — n „■ — 1 / ^ 

„7/5 r 

= _-y |V$r + Tr(-iA7). (30) 



From (jH)) we get that 



where 

™..„„^//va.)v..(.)i-r'*v(.).vfe)d.d,. (32) 

(Since the Coulomb energy is an unbounded operator one may worry about the convergence 
of the expansion in (jHT|) . This problem is easily circumvented by introducing a convergence 
factor into Ixl""*^, e.g., |x|~"'^(l — exp(— t|x|)). The expectation on the left of (j3ip converges 
as t ^ oo by the Monotone Convergence Theorem, since for fixed values of the charges each 
term is monotone in t. We may do all calculations and estimates for finite t and at the end 
let t -^ cxD. We will here ignore this slight complication.) 
Using the notation of Section |21 we have 

(V'/3e', 7lV'ae) = ^ l -X^ ^°^' ^l(^/3e'' ^«e) = 2"!^^'^°'^ ^^^^ 
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and thus from (j24j) 



«e«^e'V'"^e) = ^(00, V'a) (00, ^/j) (^/., 0o) (V^^, 0o) 

-n— ( ^»I3 ^ _"^2 (^M^ 0o) ii'u, 0o) + ^M'^ ^ _^^2 ('^O' ^") ('^O' ^/^) 



+n— I 5, 



afi 



-. ^(00, 'ipp){'ipu, 0o) + S/Bu- r2-(0O, V^a)(V^M' 00 



/3 



^2 -,2 

+ 7'^/5m1 T2"(0O' V'a)(V'i', 0o) + Jf^ai.^ Ti'^'^O, ^/3)(^/., 0o) 

4 i - A^ 4 i - A^ 



^au^lSfM ^a ^P I ^a^iS|3u A^ A^ Sa/S^^u Aq, A^ 

4 1 - A2 1 - Al + 4 1 - A2 1 - Al + 4 1 - A^ 1 - A^ 



We therefore arrive at 

(oo \ oo / \2 \ 

*-© E 7-^^n =E^"-^(^-<^o)(^.,0o)4t^-T^ 



X2 

where we have used that 0o and V'a, « = 1, • • • are real. From the expression for Waafiu we 
see that we may write this as 

' oo \ 

^ N=0 1<i<j<N ' * •'I / 

where /C is the operator on L^(]R^) with integral kernel 

IC{x,y) = (l)oix)\x-y\'^(l)o{y). (36) 

Putting together (I^U)) and (jH3)) we arrive at 

7/5 p 

(i7(2)) = ^y |V$r + Tr(-iA7) +nTr (/C (7 - ^/^(^Tiy) ) . (37) 

Our next goal is to construct the operator 7. Here we shall use the method of coherent 
states symbols. Let x{x) = 7r^'^/^exp(— x^) such that /x(a;)^dx = 1. Let < £ be a 
parameter which we shall specify below as a function of n such that n^^/^ <^ i <^ n^^/^. 
Denote Xiix) = £~^/^x(x/£) and let 

6u,p{x) = exp{ipx)Xe{x - u). (38) 

We then define 7 to be the operator 

7 = (27r)-3 If f{u,p)\eu,p){e^,p\ dn dp (39) 
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where 

We see that f{u,p) > and hence 7 is a positive semi-definite operator and since f{u,p) = 
f{u, —p) all eigenfunctions of 7 may be chosen real. That this is an appropriate choice for 
the function / will be seen at the end of our calculation (see (jlSj) ). Moreover, 

Tr7 = {27r)-' JJfiu,p)du dp = 7T''/^(^^Y' j^Muf/'duj ^gip) dp 



Thus 7 is a trace class operator. Hence we have all the requirements needed in order for 7 
to define a state \l/„. Moreover, we see from (|^ that for large n 

(AT) =n + C(ri3/5). (42) 

We turn now to the calculation of the expectation of the kinetic energy. 

Tr(-A7) = {2n)-' jj j\V9^J^f{u,p) du dp 

= (27r)-3 ffp^f{u,p) du dp+ (27r)-3 f {VxYe"^ f f f{u,p) du dp 

< (27r)-3 ffp''f{u,p) du dp + C(n2/5£)-'n^/^ 

^23/V-7/V/5 f $(^)5/2d„ f p^g(p)dp + C{n^/H)-^n''/\ (43) 



where in the second to last inequality we have used the definition (j^Uj) of 00- 

The next step in calculating the energy expectation in the state \l/„ is to calculate (or 
rather estimate) Tr(/C(A/7(7-|-l) ~7))- I^ order to do this we shall use the operator version 
of the Berezin-Lieb inequality given in ()76|) in Theorem lA.ll in Appendix ^ We will use 
it for the operator concave function ^(t) = \/t{t + 1) — t (see the discussion at the end of 
Appendix IX|) and the map u \-^ \u) being {u,p) ^— *> \Ou,p)- We have 

{27i)-' f\9u,p){9u,p\dudp = I. 
Since /C is a positive operator we conclude from Theorem lA. II that 

Tr(/C(V7(7 + l)-7)) 
> (27r)-3 J J [^/fiu,p)ifiu,p) + 1) - /(«,p)) {9^,p\}C\e^,p) du dp. (44) 
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Since |x — ?/| ^ is a positive definite kernel we have for < 6' 
{Ou,pmOu,p) = / / e'P'^Xiix - u)(t)o{x)\x - y\'^e-'PyXi{y - u)4>o{y) da; dy 

> (1 - C6')M^? ffe'P^Xiix - u)\x - y\-'e-'PyXi{y - u) dx dy 
-C5''\n"HYn-^"> 

> Muf 11 e'''^Xi{x)\x - y\-\-'^^Xt{y) dx dy - Cb\n^lHfn-^^'' 
-C5'-\n"Hfn-^"> 

> Mnf I jM)t^, dg - C{n'/Hfn-^/\ (45) 

J \P-(1\ 

wliere jn^q) = (27r)~^|X£(g)p = £'^7r~^e"^^ '^ (witli tlie convention f{p) = /e*^^/(a;) dx for 
tlie Fourier transform). In the last inequality we have chosen 5' = {n'^^^i)n^^^^ and in the 
first inequality we have used that \4>o{x) — 0o(^)| ^ Cn^^^\x — u\ and hence 



Xi{x - u)\4>o{x) - 0o(m) \\x-y\ \i{y- u)\My) - 4>o{u) \ dx dy < C{n^/HYn ^/\ 
We have that f ji{q) dq = 1. We will use the estimate 

/* II /* II 

/ \q\ I \q\ 

< l^r^ / Jiil)] r dq+ IpI"^ / je{q)-. ^ dg 

J \p-q\ J \p-q\ 

< sup {Uq)\q\'^') (^\P\-' J \q\-'^'\p - q]'' dq + \p\-' (^j \q\-''''\p - q\-' d<^^ 

< Cbr^/^up(j,(g)|gr/^). (46) 

For our explicit choice of j^ we get \\p\~'^ — ji* bl^^l < t^^^'^\p\^^/'^. From (HH), (jl^j) and 
estimate PT|) we find that 

Tr(/C(v/7(7 + l)-7)) 
> 2(27r)-2 jj {yf\u,p){f\u,p) + l) - f\u,p)) Mufji * br' du dp 



> 2-i/V-^/V/^ jj iV9{p){9{p) + l) - 9ip)Muf'\pr' du dp 

- C{n^/'i)-'/^n^/' - C{n'/H)^n'/', (47) 

where we have also used that // ( ^/fiu,p){f{u,p) + 1) — fiu,p)] du dp < Crr'^^ (as in 

m)- 
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If we now insert the above estimate and (j43j) into (j37j) we arrive at 
(i7(2))<^7/5/'i f |v$(^)|2dn 

^2-l/V-7/4 f $(^)5/2^^ f p^g(p)-\p\-^Ug(p)(g(p) + l)-g(p)) dp 

The function g in ()40|) was chosen precisely so as to optimize the above expression. If we 
insert the expression for g it is easily seen that the term in the large parenthesis above is 



JR3 J]R3 



If we choose $ to be an exact minimizer then this expression is —A (recall that A and Iq 
were defined in Theorem ll.2p . From the estimate in (|48p we see that if we choose £ as a 
function of n such that in^^^ = v?!'^'^ then 

(7f(2)) < _An^/5(l _ Cn-^/35)^ (4g) 

Because of the estimate (jl^ this means that we have found a state satisfying (j^3j) . 

We could instead have chosen $ to be a smooth compactly supported approximate min- 
imizer to the variational problem (@)). We would then for any e > have proved that 
lim„^oo '^~ (-f^ ) < —A + e, which of course implies fl25p . 

3.1 An upper bound for fixed particle number 

In this section we shall prove the upper bound in Theorem II .21 on the energy E^'^\N) corre- 
sponding to a fixed particle number A^. 

Let "^e^n for n,e > denote the state constructed in the previous section, but with the 
function g in (j4Uj) replaced by the function g^, which is equal to g for \p\ > e and is zero 
otherwise. We will again denote the expectation of any operator A in the state \E'e,n, by (A). 
It then follows from the construction in the previous section that 

lim n-^/^(/7(2)) < _A^^ (50) 

n— >oo 

where A^ — > A as e — >■ 0. 

Let \l/e^ denote the projection of the state '^e,n onto the subspace corresponding to 
particle number m = 0,1, . . .. We then have 

m=0 \ \e=it «=1 
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Hence from ^ and (jMl) 

oo 

{Af^) - (7V)2 = J] 5^ «e«a,e«L'«a,e') " ("La^.e) «e'aa,e') =n + 2Tr7^(7^ + l), 
a=l e,e'=± 

where 7^ is given as in ()39|) . but with / replaced by /e, which is expressed in terms of Qg 
instead of g. Thus using (I75|) in Theorem lA.il (or (j76p for that matter) in the convex case, 
we see that 

(Af^) - (Af)^ <n + 2(27r)-3 // Uu,p){Uu,p) + 1) du dp < n + C,rv'l\ 



Here C^ > is a constant depending on e and such that C^ ^ 00 as e ^ 0. It is at this point 
that it is necessary to replace g with ^f^, since otherwise the above integral is not convergent. 
For any M > we have 

m-{M)>M m=0 

< M-3/5(Ar2)ViO((Ar-(Ar))2)3Ao 

Given a positive integer A^, we choose n = N — CqN^^^ . Then if Co > is chosen 
appropriately we have according to (j^ and (jl^ that the expected particle number satisfies 



N - CiN^^^ <{Af)<N- C2A^^/^ 

for some Ci, C2 > 0. 

Since M \-^ E{M) is a non-increasing and non-positive function (adding particles will 

always lower the energy, since one may construct a trial state with the extra particles placed 

arbitrarily far away from the original particles) we have that 

00 
i^(^)(Ar)<5^E(2)(m)||vI/(-)f < $:i?(^)(m)||v^(-)f- Y. ^^'HHII^Sf 

where we have used the lower bound E^'^\m) > —Cm7^^ (see jS] orJT]) and the estimate 
^^. Thus we finally get the upper bound in Theorem 11.21 



limsupA^-^/^E(2)(A^) < limlimsupn~^/^((i7(2)) + C,A^^/^"^/^°) = -A, 

jv— ►oo e— >U n— >oo 

according to (j^ . 
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4 The one-component charged Bose gas 

Since the thermodynamic ground state energy e(p) of the one-component charged Bose gas 
may be calculated by minimizing over all particle numbers we may again consider the grand 
canonical ensemble. Thus we are looking for an upper bound to the ground state energy of 
the Hamiltonian H^^' = ©^=o-^iv acting on the Bosonic Fock space JF(L^(A)). 

To construct a grand canonical trial function we begin by choosing a real normalized 
function 0o ^ -^^(^)- Let i] G Cq(0,L) be a non- negative function compactly supported in 
(0, L) and such that L ri{tY dt = 1. Moreover, assume that ri{t) is a constant for t G [r, L— r] 
for some < r < L/A to be chosen below. We will write this constant as {p/nY'^, for some 
n > 0. In fact, we shall choose r independently of L (for large L). We also assume that 
v{t) ^ {p/ny^^. We then define 

(j)oix, y, z) = r]{x)r]{y)r]{z). (52) 



Thus 00 is equal to a constant \fpjn on the cube [r, L — rf and < 4>q{x) < ^/p/n for all 
x G A. Since rj is normalized so is (pQ and p{L — 2r)^ "^ n < pL^. Thus the constant n is 
almost the number of particles required to have a neutral system. We have 

\r]{t)\<CL-^/^ and \(f)o{x)\ < C L^^^^ (53) 

and we may assume that the derivatives satisfy 

W{t)\<Cr-^L-^/^ and hence \VMx)\ < Cr-^L-^^^. (54) 

In particular, we have 

/|V0o(a;)pdx<C(rL)-^ (55) 

Ja 

Observe that we also have that 

{n(f)o{xY — p)\x — y\^^{n(f)o{yy — p) dx dy < Cp^L^r^. (56) 

We choose our grand canonical trial function \I'„ as in (j22|) . The condensate vector is 

= 2^000 (57) 

where the parameter Zq > will be chosen below. The operator 7i = 7 (we omit the subscript 
1 because we shall use a subscript e below with a different meaning) will be chosen to be 
a positive semi-definite trace class operator with real eigenfunctions. The eigenfunctions 
(corresponding to non-zero eigenvalues) should satisfy Dirichlet boundary conditions on the 
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boundary of A. Let ipa, a = 1, ... be an orthonormal basis of real eigenfunctions for 7. We 
use the notation a* = a*{ipa)- 

As usual we denote the expectation of an operator A in the state \l/n by {A) . As in (jHUj) 
we see from (0) and (j^ 

^ 00 N \ 2 r 

*- E -i^^*" = f / '^'^°'' + Tr(-lA7) < CzlirL)-' + Tr(-lA7), (58) 

^ N=0 i=l ) -^ 

where in the last inequality we have used (fS^ . We likewise get 

^ 00 Af \ c c 

*- Y. ^(a^O^n = / V[x)<\>{xf dx+ V{x)p^{x) dx 

^^ — n ^ — i / ^ ^ 



N=0 i=l 

^2a („,\2 



= p II 'oMyf + P-^^y^ dx d,, (59) 

JJax/v \^ ~ y\ 

where p^(x) = 7(x, x) is the density of the operator 7. 

From (jHl) we have (as in EI} with Waisufj. given exactly as in (jH^ 

N=0 l<i<j<N / a,l3,tj,,u=l 

We then obtain from ()24|) that 

^"'0 Zl l^i-^J-r'^") =y // 0o(a;)'|x-i/rVo(l/)'dxdi/ 
+ZqTt [JC (^7- a/7(7 + 1)JJ +4 M^fl^-yl'^Pyi^) da;dy 

JJaxa \^ " y\ JJaxa f — y| 

+ 1// P7(3^)k-Z/rV7(l/) dxdy, (60) 

J JAxA 

where the operator /C is given as in (jHBj) . Putting together (j3Hj) . (j3^ . and (jUUj) we arrive at 



AxA 



+Tr (-iA7) + Zo'Tr (/C (7 - ^7(7 + !))) • (61) 

We now choose 

^ = ^' = <'''"I*((4w)l'''><*^l* (62) 
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where the function gs{p) = for \p\ < e and ge{p) = g{p) for |p| > £ where g is defined in 
p]|l . and 

9p{x) = \/np~^ exp{ipx)(f)o{x). (63) 

Recall that np'^(f)Q{x)'^ < 1 and is equal to 1 on most of A. 

We see that the map p \-^ \6p) satisfies the requirements of the map uj ^—>- |u;) in Theo- 
rem ^3 with measure djjiiuj) = (2-7?) ~^ dp. 

That 7e satisfies the necessary requirements follows as before. It is clear that the eigen- 
functions of 7^ with non-zero eigenvalues have compact support in (0, L)^. 

We calculate the density of % 



P.(x) - i^^r' l9^[j^. ]K(-y\^p 



(27r)-V-Vo(x)^ / gs { j^^, ] dp 



= np-^/^2-='/V-9/Vo(a;)' / geip) dp. (64) 

We finally choose Zq > 

zl = n(l- 2-'/'p-'/\~'/' J gM dp) (65) 

(for p large enough). Then 

Zo4>o{xf + p^,{x) = n(j)Q{xf. 

It follows from (|56|) and the fact that (t)Q{xY < p/n that 

(p - P7.(^) - zlu^f) \x - y\'' (p - PiAy) - zlUyf) dx dy < Cp'LW (66) 

AxA 

To estimate the second term in (pT|) we will use Hardy's inequality / |V'u(x)p dx > 
\ J \^ dx as follows 



X — y\ \j J / \j J \x — y\ 

1/2 / ^ ^ \ 1/2 

< 2 ' ' ' 



// be{x,y)\^ dx dy] I // \V^%{x,y)\^ dx dy 
2(Tr7,^)^/^(Tr(-A7,^))^/^ 



Since a; i— > x^ is operator convex we may estimate these terms using the Berezin-Lieb in- 
equality (fTBj) in the convex case, but we may alternatively simply use the norm bound 
||7e|| < Ce'"^. Hence 

\x-y\ 



dx dy < Ce^^ (j P,A^)\ (Tr(-A7,))^/2 < Ce-^pL\ (67) 
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where we have used (j64j) . n < pL^ and the, fact which we shall prove below in (j68|) . that 
Tr(— A7£) < Cp^/^Li^ (recall that we will choose r independently of L). The third term 
in ()61|) which compared to the second term has 7^ replaced by a/tTOtT+I) is estimated in 
exactly the same way and with the same bound as the second term. 

We are now left with calculating the last two terms in (jF)T|l . For the kinetic energy of 7^ 
we have as in (jiHj) 

Tr(-A7.) < (2.)-^/^^,.(^)(/ + /|VM^)rd.)d. 

< 23/47r-^/V/4L3 / pV(p) dp + Cp3/4i:3(rL)-i/2, (68) 

where we have used (|55|) and n < pL^ . 

For the last term in (jfilj) we again, as in ()44|1 . appeal to the operator version ()7fi|l of the 
Berezin-Lieb inequalities. We arrive at 

Tr(/C (7,-^7.(7. + !))) 

< (27r)-3 / (/,(p) - V/e(p)(/e(p) + l)) {e,\lC\\e,) dp, (69) 

where /e(p) = Qe {pi^T^p)^^^^)- We have as in (jIH|) 

(^p|X:|^p)=47^J*br^ (70) 

where J(p) = (27r)~^np~^|0o(p)P- The special form (j^ implies that 

J{Pl,P2,P3) =j{Pl)j{P2)j{P3), 

where j{t) = {27i-)~-^n^/^ p~-^^^\ri'^{T)\'^ . Since J J{t) dr = n^/^p^^^^ j rj{tY dt, /?7^ = 1, and 
< ''7(t) < n^^^'^p^^'^ and equal to this constant on [r, L — r] we have that 1 — 2r / L < 
J J{t) dr < 1. This implies in particular that 

(1 - 2r/Lf < I J{p) dp < 1. (71) 

By (|53p and (j54|) and the support property of f]' we we have |?7^(t)| < |r|^^ J |(^^)'(t)| dt < 
C{\r\L)-\ Thus j(^) < CL{\T\Ly\ Hence 

/ J(g) dg < 3 / jir) dr < CL-V2. (72) 

^|g|>L-i/2 J|r|>(3L)-V2 

For IpI > £:(87rp)-'^/'^ and |g| < L~^^'^ we have \p — q\ < (l + C p^^/^e^^L^^^'^)\p\ and hence 
from (HH) and ((721) 

J*\p\-^ > {l + Cp-^/^e-'L-'/Y^\p\-^ [ J{q)dq 

J\q\<L-^/^ 

> (1 + Cp-^'^e-^L-^'Wil - 2rL-i)=^ - CL-^'^)\p\-^ 

> (l-C(p-i/ViL-V2 + rL-i + L-V2))|p|-2. (73) 
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Inserting this into (j7(jp and then into (j69p we arrive at 
Tr(/C(7.-V7e(7e + 1))) 

< 2-^"p^'''K-"' j{g,{p) - ^ge{p){ge{p) + mp\-' dp 

+C{e-^L-^'^ + p^/Vl-i + pi/'^L-i/^). (74) 

If we now insert the above estimate, (|H^. (jUUI) . (jUH|) . (jUTj) . and the same estimate for 7^ 
replaced by v^7e(7e + 1)) into (jnij we see that 

hmSUpL-3(/7«) < p5/42-l/4^-7/4 /■ |p|2^,(p) +^,(p)|p|-2 _ ^ g,{p){g,{p) + l)|p|-2 dp 
L— >oo J 

+Cp(l + pr2 + £-2). 

Here we may actually let r — > (which really means that we could have chosen r as a 
negative power of L). If we recall the behavior of g{p) for small |p| from (PUJ) we find that 
the error in replacing (7^ by g is of order p^l^e. Thus by choosing e = p^^/^^ we obtain the 
final result 

e(p) < \imsnp L'^H^^^) < -hp'''\l - Cp-^'^^). 

L— >oo 

A The Berezin-Lieb inequality 

In this appendix we shall prove variants of the Berezin-Lieb inequalities j21 E] ■ 

THEOREM A.l (Berezin-Lieb inequalities). Let H be a Hilhert space and Vt a measure 
space with a (positive) measure p such that there exists a map 

VL3 oj ^^ \uj) ETi, 

satisfying J \uj){uj\dp{uj) < I as operators. Assume ^ : M_|_ U {0} ^ M. is a concave function 
with C,{0) > 0. Then for any non-negative function f onQ satisfying J f{u!){uj\uj)dfi{uj) < 00 
we have the Berezin-Lieb inequality 

Trn U (J f{uj)\u){co\dfi{cu)\\ > j ^{f{uj)){uj\uj)dp{u:). (75) 

If moreover ^ is operator concave (still satisfying ^{0) > 0) the inequality holds as an operator 
inequality 

e ( [ f{u;)\u;){uj\dfi{uj)] > /e(/M)k)(^M/iM. (76) 
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Proof. We first note that J f {uj)\uj) {uj\dfi{uj) is a positive semi-definite trace class operator. 
Let Ui,U2, ■ ■ ■ be an orthonormal basis of eigenvectors for tliis operator. Tlien 



Tr„ U / fiLu)\cu){iu\df^iu) = yj[ fiuj)\{u\u,)\'dfiiuj) 






oo 



i=l ^ 
xd ( [ \{u\uM'dfiiLu) ] I fiLu)\{iu\u,)\-'df^icu) 1 , 



wliere we liave used tliat J \{uj\ui)\'^dfi{uj) < 1 and tliat since ^ is concave witli .^(0) > we 
liave ^(at) > aC,{t) for all t > and < a < 1. If we now use Jensen's inequality we arrive 
at 



We turn to the case when ^ is operator concave. Define the operator U : 7i -^ L'^{Q, dfi) 
by (t/0)(cu) = (tu|0). Then 

U*h= f h{uj)\u)dfi{uj). 



Thus if B is the multiplication operator on L'^{Q, dfi) given by Bh{uj) = f{uj)h{uj) we have 

U*BU= I f{uj)\uj){uo\d^i{uo). 



In particular, we have the operator inequalities < U*U < I. Using that (1 — UU*Y^'^U = 
f/(l — U*Uy^'^ it is straightforward to check that the following operators on 7i © L^(r2, dfi) 
(written in matrix notation) are unitary 

{i-u*uy/'^ -u* \ ({i-u^uy^ u* 

U=\ ' ^= 

u {i-uu*Y'^ I \ u -{i-uu*f/^ 

Moreover we have that 

2lofiy 2I05/ lo (1- uu*f/^B{i ~ uu*Y/^ 
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Since ^ is operator concave and U and V are unitary we find that 

i{U*BU) 

^((1 - f/f/*)i/25(i _ f/f/*)l/2) 

i*/o o\ i,/o , 

U*i{B)U 

{i~uu*Yi^i{B){i-uu*Yi^ 

In particular, this gives E,{U*BU) > U*^{B)U, which is precisely the operator Berezin-Lieb 
inequality (f7n|l . D 

In order to determine whether a given function is operator concave we may use Nevan- 
linna's Theorem (see [3] Theorems V.4.11 and V.4.14 and equation (V.49)). According 
to this a real function ^ defined on the positive real axis with an analytic extension to 
C\{xgM|x<0}, which maps the upper half plane into itself has a representation of the 
form 

where /3 > and where i/ is a positive measure satisfying J^ iTj^d^W < oo. Since 
t I— > — (t + A)^^ is operator concave the same is true for functions with the above integral 
representation. 

As a special case we see that the function ^(t) = ^/t{t + 1), which is analytic away from 
the segment [—1,0] is operator concave. 
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